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Classical-like behavior in quantum walks with inhomogeneous, time-dependent coin 
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Although quantum walks exhibit peculiar properties that distinguish them from random walks, 
classical behavior can be recovered in the asymptotic limit by destroying the coherence of the pure 
state associated to the quantum system. Here I show that this is not the only way: I introduce 
a quantum walk driven by an inhomogeneous, time-dependent coin operator, which mimics the 
statistical properties of a random walk at all time scales. The quantum particle undergoes unitary 
evolution and, in fact, the high correlation evidenced by the components of the wave function can 
be used to revert the outcome of an accidental measurement of its chirality. 
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I. INTRODUCTION 

Quantum walks (QWs) N were originally termed 
“quantum random walks” as they were thought 

as the quantum-mechanical version of the discrete ran¬ 
dom walk (RW) in one dimension: the Markov process in 
which a particle changes its position at each clock tick by 
jumping to one of the two nearest sites depending on the 
random outcome of a coin toss. This source of random¬ 
ness could be seen as superfluous in the quantum world, 
where the location of a particle is a probabilistic magni¬ 
tude, governed by its wave function. Therefore, in the 
design of these “quantum random walks”, the coin toss 
was replaced by some (unitary) operator that affects the 
state of a quantum binary property of the system, e.g., 
the spin or the chirality, and the wave function is shifted 
according to the value of this qubit. 

Consequently, beyond the intrinsic uncertainty of the 
quantum phenomena, “quantum random walks” are not 
random at all —and thus this term is now deprecated. 
The most prominent sign of this deterministic nature 
of QWs is the ballistic behavior they can show the 
ability to connect any two sites after a lapse of time 
that is proportional to the distance between these sites, 
even if the walk is undirected. This fact comes in con¬ 
flict with the diffusive nature of unbiased RWs which, 
to perform the same operation, need a time period that 
grows quadratically with the separation of the sites. This 
speed-up readily caught the attention of the scientific 
community, albeit there are other properties that distin¬ 
guish QWs from RWs Q- In spite of those differences, 
QWs are indeed the quantum analogues of RWs, and 
therefore they may experience a change from ballistic to 
diffusive motion when the quantum coherence of the state 
is affected by multiple reasons [sl-fT^. In fact, it has been 
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proven that, under mild conditions, the introduction of 
temporal or spatial fluctuations in the properties of the 
coin operator acting upon the QW leads to classical be¬ 
havior in the asymptotic limit: the standard deviation of 
the position of the quantum walker grows with the square 
root of the elapsed time, and the corresponding rescaled 
distribution converges to a Gaussian [T^ . 

Soon after the birth of the very concept of quantum 
computers [l^, i.e., computers whose operation can¬ 
not be understood without the laws of quantum me¬ 
chanics Il6|. t he first genuine quantum algorithms ap¬ 
peared [l7Hl9l| . algorithms that were more efficient than 
their classic counterparts. And since many of those clas¬ 
sical algorithms use RWs as building blocks, it is not sur¬ 
prising that the ballistic transport of QWs was seen as 
the key feature in the design of faster algorithms [20l - [2^ . 
But QWs can play an even more important role in quan¬ 
tum computation, as they may be regarded as universal 
computational primitives [1^ 241, i.e., they can be used 
to implement all the logic gates that a universal quantum 
computing machine needs to work 

While it is still an open issue whether universal quan¬ 
tum computers can efficiently simulate arbitrary physical 
systems |26| . in this paper I will give an affirmative an¬ 
swer to a related but much more limited question: Can 
a quantum walk be used to simulate the behavior of a 
classical system whose evolution follows a random walk? 

Specifically, I will look for a QW that shows exactly 
the same probabilistic properties of a RW at all time 
scales. This objective must not be reached as a result 
of the introduction of exogenous disturbances that can 
induce decoherence in the pure quantum state: The pur¬ 
pose is to obtain the classical distribution (the binomial 
distribution) by means of reversible unitary evolution at 
every time step. Therefore, this QW with classical-like 
attributes could replace the corresponding RW in the 
simulation of the dynamics of the desired classical sys¬ 
tem on a quantum computer. 

With this aim, I consider here a discrete-time QW 
on the line endowed with an inhomogeneous, time- 
dependent coin operator. Extensions of this kind were 
considered in the past: one can find in the literature ex- 
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amples of QWs driven by inhomoKeneous, site-dependent 
coins [27l - K^ . time-dependent coins [32l - K^ . or history- 
dependent coins [13, 

The paper is organized as follows. Section |TI] reviews 
the formalism used in the construction of the discrete¬ 
time quantum walk on the line with a time-dependent 
coin operator. Section mil shows how one can devise a 
QW that behaves like an unbiased RW. I extend the 
framework to encompass general RWs in Sec. IIVI Sec¬ 
tion El explores the possibility of reverting the conse¬ 
quences of a measurement of the chirality, and quanti¬ 
fies the entanglement between the chirality and the posi¬ 
tion of the particle by means of von Neumann’s entropy. 
When this correlation is destroyed. Sec. m the walker 
becomes vulnerable to decoherence. The paper ends with 
Sec. IVIIl where conclusions are drawn. 

II. QW WITH AN INHOMOGENEOUS, 
TIME-DEPENDENT COIN OPERATOR 

I begin the discussion by introducing the foundations 
of the inhomogeneous, time-dependent quantum walk on 
the line. I denote \yy TiP the Hilbert space of discrete 
particle positions in one dimension, spanned by the ba¬ 
sis {|n) : n G Z}, and by 'W the Hilbert space of coin 
states, spanned by the basis {|+),|—)}. The discrete¬ 
time, discrete-space quantum walk on the Hilbert space 
= the j;esult of the action^f the evolution 

operator Tt, Tt = SUt, where the coin Ut is an inhomo¬ 
geneous, time-dependent, real-valued unitary operator: 

CXD 

[cos6»„_t|-b)(-b|-f sin6»„^t|-h)(-| 

n——oo 

+ sin6»„_t|-)(-b| - cos6»„_t|-)(-|] 0 |n)(n|, (1) 

with 0 < ( < TT, and S is the shift operator that moves 

the walker depending on the respective coin state: 

S'li) 0 |n) = |±) 0 |n ± 1). (2) 

As the time increases in discrete steps, one chooses the 
time units so that the time variable t is just an integer 
index, and the state of the system at a later time, 
is recovered by applying Tt to the present state 

\i;)t+i=ft\^)t. (3) 

Equation m induces the following set of recursive equa¬ 
tions: 

+ = cosO^ t -tsin6»„_t 

(4) 

- cose^t 

(5) 

on the wave-function components, the projec¬ 

tions of the state of the walker into the elements of the 


basis of the Hilbert space: 

ip+{n,t) = {n\<Si {+\ip)t, ( 6 ) 

= {n\<Si (7) 

The evolution of the system is fully determined once 
\ip)o = |V’)t=o is set. Since the final aim is to simulate a 
RW, we must consider that the particle is initially located 
at the origin. When the coin operator is homogenous and 
time-independent, it is well known that the chirality of 
such localized state affects the ulterior behavior of the 
system [Ulil. In our case, as we will see later on, this 
choice is not so delicate. Thus, for the sake of simplic¬ 
ity, I assume that there is no preferred direction in the 
chirality: 

IV’)o = ^(|+) + |-))®|0), (8) 

that is '!/’±(0,0) = l/-\/2- Note that a real-valued state at 
time t = 0 precludes the possibility of having a complex¬ 
valued wave function at a later time, cf. Eqs. and ©■ 
We want to connect the evolution of our quantum sys¬ 
tem with the statistical properties of a random walker. 
This connection must be done through the pairing of the 
probability mass function (PMF) of the two processes: 
Let us call p{n, t) the likelihood of finding the particle 
in a particular position n at a given time t. In the case 
of the quantum walker this probability depends on the 
wave-function components, 

p{n,t) = |V’+(n,t)|^ -b . (9) 

Our goal is to get that p{n, t) equals the PMF of a random 
walk, the binomial distribution: 

pM = (1 -P)^ , (10) 

for n G {—t, —t -b 2, • • • , t — 2, t}. Function p{n, t) deter¬ 
mines the different moments of the stochastic process, 

t 

n— — t 

among which are worth to be highlighted the expectation 
value of the walker position, (A)^, and its uncertainty 
AAj, magnitudes that should amount to 

{X)t = {2p-l)t, (II) 

AA, = \/(A2),-(A)2 = 2^p{l-p)t, (12) 

if the classical expression is valid. Note that, since 
the probability is conserved, the change in the expecta¬ 
tion value of the position at consecutive instants reads 

t 

(A),+i = (A),+ X •^(^.i)> (13) 

n——t 
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where J{n^ t), 

J{n,t)= cos 26»„ t [ipl {n,t)- ip^_ {n,t)] 

+ 2sm2e^^tj}_^_{n,t)1p_{n,t), ( 14 ) 

is the net flux of probability entering or leaving site n, 
and its explicit expression stems from Eqs. o and ©■ 
As we will see in the next section, Eqs. m and m 
pave the way for achieving our purpose. 

III. UNBIASED WALK 


Note that for n 7 ^ 0, ^ In fnct 

tp^{n,t) = - 2,t), see Fig. [U a property whose 

implications I discuss below. Once one has the explicit 
expression for the components of the wave function, the 
coin weights read 

sta9„,, = l(v/l+|^+V'^^). (26) 


Our task is therefore to deduce a functional form for 
cos 9^ ( that can be accommodated in Eqs. and m 
and ultimately lead to the desired PMF, Eq. (HOD . In 
order to grasp the appropriate procedure, I will consider 
the unbiased version of the RW in the first place, 

1 t\ 


for n G {—t, —f + 2, • • • ,t — 2,t}. This results in a great 
simplification since in this case the expectation value of 
the position is null , (A)^ = 0, for any time value. This 
property is preserved by Eq. 0 if J{n,t) = 0, a suffi¬ 
cient condition. The absence of probability flux can be 
readily achieved, see Eq. (O, if 


cos 26»„ ^ = 


sin 29^ t = 


2 V)_^(n,t)V)_(n,f) 

p{n,t) 

p{n,t) 


( 16 ) 

( 17 ) 


that is. 


cos6i„t = 


sinfi'n,* = 


1 V’+(rt,f) - V>_(n,t) 
\/2 ^/p{n,t) 

1 ip+{n,t) + 

\f2 p(n,t) 


( 18 ) 

( 19 ) 


It is easy to check that Eqs. (na and (HU) represent valid 
trigonometric expressions. Now, one can introduce these 
formulas in Eqs. and ® and obtain: 


il^^{n + l,t + l) = 'ip_in- l,t-k 1) = y ^^”’^\ (20) 
leading to 

for n G {—t + 2,—t + 4, - ■ ■ ,t — 4,t — 2}, and 

V’+(t,t) = V'-(-^,t) = , (23) 

='ip_{t,t) =0. (24) 
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Figure 1. (Color online) The two components of the wave 
function at t = 31. The red dots correspond to ip+{n,t) 
whereas the blue diamonds mark the values of ip_{n,t). 


IV. BIASED WALK 

All these results can be easily modified to encompass 
the generic case: we simply need to replace the factor 2 “* 
in Eqs. m and (I22p by the proper combination of powers 
of p and (1 — p). Moreover, conditions (pi?)) and (IMll 
should be mapped into 

'>p+{-t,t) = = 0 , 

= {1-py , 


what suggests the choice 


^+in,t) = (l-p) ^ , (27) 

V'-KO = (1 1 (28) 

for n G {—t -I- 2, —f -I- 4, • • • ,t — 4,t — 2}, see Fig. [21 In 
other words, Eq. (1201) now splits into 

'fp+{n,t) = 


(29) 

(30) 
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Finally, we have to use recursive Eqs. and m to iso¬ 
late cos0„ j and sin0„ 



which satisfy all the desired constraints. Note how ex¬ 
pressions (ED and (1321) are ill defined for n = t = 0: in 
fact, Eqs. m and m evidenced this same issue. To be 
consequent with the previous setup and, in particular, 
with Eq. (IHl), the right option is the most obvious, i.e.. 



but since our coin operator is time dependent, one could 
modify 0g q and at will, as long as one has 


li’)! = VpI+) ® |i) + \/i -p|-) ® I - 1), (35) 

unchanged. This invariance is just one of the many possi¬ 
ble transformations that preserves the functional form of 
p{n, t) [4l|, but this will be the subject of future research. 



Figure 2. (Color online) The two components of the wave 
function at t = 100, for p = 0.75. The red dots denote ip+ iji, t) 
whereas the blue diamonds designate the values of 


V. REDUNDANCY AND COHERENCE 

Consider now the following prominent consequence of 
Eqs. (1291) and (l30l) . On the one hand, we have a high 
degree of redundancy, with almost the same information 
stored in each component of the wave function. On the 
other hand, this information is the PMF of the system 
one time step before. All together implies that one can 
undo the consequences of an accidental measurement of 


the chirality at time t, by means of unitary transforma¬ 
tions. In particular, if |'0)t —>■ one has 

*) = (1 - P) ~ > (36) 

^l(n,t)=0, (37) 

and the recovery procedure is 

\fj)t=LSV+\^+)t, (38) 

where V~^, 

[Vp|+)(+| + v^I+)(-| 

+ V^|-)(+|-V^|-)(-|]0/P, (39) 

is a homogeneous coin operator, and L, 

CO 

L = r® ^ |n—l)(n|, (40) 

n— — oo 

represents a systematic shift to the left. Thus, the joint 
operation of LS displaces the negative component of the 
wave function two sites to the left, whereas the positive 
component remains in place. On the contrary, if one has 
obtained the unitary operation is 

H)t=RSV-\^-)t, (41) 

with 

[-yr^i+)(+i + v^i+)(-i 

+ v^|-)(+l + V^|-)(-|]o?", (42) 

and 

OO 

^ |n-|-l)(n|. (43) 

n— — OQ 

In the last expressions P and P denoted the identity 
operator of the corresponding Hilbert space. 

The procedure just described can revert the system to 
the previous state provided that the outcome of the acci¬ 
dental measurement of the chirality is known. Otherwise, 
the quantum walker will suffer decoherence since there is 
some probability that one chooses the erroneous unitary 
transformation, i.e., that one applies Eq. (1381) when the 
right choice is ED, and vice versa. 

Consider, for instance, that the fortuitous measure¬ 
ment occurs after the first time step yielding 

l^“)i = I-) O I - 1), 

and that we perform the wrong unitary transformation: 

\cj>)i=LSV+\^-)i 

= \/l -p|+) o I - 1) - Vp|-) ® I - 3)- (44) 
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The new PMF is clearly different from the one we had 
but in a way that is difficult to quantify. Standard meth¬ 
ods of information theory, as the Kullback-Leibler diver¬ 
gence [3, require that the null sets of both probabil¬ 
ity measures are equal, which is not the case here, cf. 
Eqs. (n?5l) and (PI) . One cannot resort to simple quan¬ 
tum arguments to establish the resemblance between \(j))t 
and I'!/')( either: Note that as 

('0+|'0“)t = 0, 

one must have necessarily 

{'4’\4>)t = 0 , 

since the recovery procedure is unitary. 

The potential relevance of an incidental and unnoticed 
measurement of the coin state will drastically depend on 
the actual correlation between the two components of the 
wave function: the higher the correlation, the lower the 
consequences. In other words, the impact will decrease 
with the actual level of entanglement between chiral and 
spatial degrees of freedom of the original quantum state. 

A way to assess this level of entanglement of the walker 
is through the entropy of entanglement [H, , a special 

instance of von Neumann’s entropy. The von Neumann 
entropy Sit) of a quantum system is defined in analogy 
of the Gibbs entropy by: 

S{t) =-tv{pt\og 2 Pt), ( 45 ) 

where pt is the density matrix operator at time t, and 
tr (•) is the trace, e.g., 

OO 

tr(Pt)= (^1 ® [(+IP‘I+) + ® l«)- 

n— — oo 

( 46 ) 

In our case, since the time evolution before the accidental 
measurement is unitary, we will have 

Pt = ( 47 ) 

and consequently S{t) = 0. However, as we are interested 
in quantifying the entanglement intensity between chiral¬ 
ity and position, one can use the reduced von Neumann 
entropy: 

5 ^(t) = -tr=(%log 2 %), ( 48 ) 

where pj is the reduced density matrix operator obtained 
when a partial trace is taken over the positions: 


( 49 ) 


with [ 45 I. li^: 

t 

^ |V^+(n,t)|^ (50) 

n——t 
t 

X] , (51) 

n——t 
t 

^ V^+(n,t)'!/'_(n,t), (52) 

n——t 

and tr'’(-) is the trace restricted to Here 

P+it)=p, (53) 

P_{t) = l-p, (54) 

for all t, see Eqs. (E3) and dsni), and 

t 

Q{t) = ^p{l-p) ^ yip{n - I, t- l)p{n + l,t- 1). 

n——t 

(55) 

Under these circumstances, the entropy of entangle¬ 
ment can be expressed in terms of (t) , the eigenvalues 
of the reduced density matrix at time t, 

S<={t) = -A+(t)log2A+(t) - Al(t)log2 AfL(t), (56) 

with, 

A^± (<) = i ± ^l-p{l-p) + Q^t). (57) 

Figure [3] shows the values of the entropy for t > 1, 
for the two examples considered in previous sections, the 
unbiased walk, p = 0.5, and the biased one, p = 0.75. In 
both instances, the entanglement is maximal for t = 1, 
when there are no off-diagonal terms in pj, and one has a 
one-to-one equivalence between the information carried 
by chirality and the spatial location —this just the case 
of the explicit example shown above, Eqs. (|35|) and (1441) . 
After that point, the magnitude of the entanglement in 
both cases converges and decreases monotonically toward 
zero: there is less information susceptible of getting lost. 

This conclusion can also be derived from the direct 
analysis of Qit) for t ^ 1. In this regime, one can ap¬ 
proximate 

p(n - 1, t - 1) ^ pin + l,t- I), 

so that 

lim Qit) = y/p(I -p), (58) 

t —¥00 

and thus 

lim S’^it) = 0. (59) 

t—^OO 

A more detailed analysis reveals that the leading term 
of the reduced entropy is of the form 

S%t) ^ ^log^it, (60) 

and therefore it does not depend on the value of p, as it 
can be seen in Fig. S) 


P+{t) = 
P-{t) = 
Qit) = 


rt = p+{t)\+)i+\+Qit)\+)i-\ 

+ Qit)\-){+\ + P.it)\-){-\ 






6 


■ . 


0.0 k 
0 



Figure 3. (Color online) Evolution of the entropy of entan¬ 
glement as a function of time for p — 0.5, green boxes, and 
p = 0.75, magenta triangles. Observe how the entropy de¬ 
creases monotonically in both cases. 


Figure 5. (Color online) Probability mass function p{n, t) at 
t = 100, for p = 0.75. Red dots were obtained from the wave 
function of the coherent process whereas the black dashed line 
corresponds to a Gaussian approximation, Eq. m- 



Figure 4. (Color online) Decay of the entropy of entanglement 
as a function of time. The green solid line corresponds to 
p = 0.5, whereas the magenta dashed line shows the behavior 
of the entanglement when p = 0.75. The black dotted line 
(~ [log 2 4t]/4t) serves as a guide for the eye. 

VI. DECOHERENCE 

From the analysis above, one can conclude that the 
present QW is particularly resistant to decoherence, but 
by no means it is immune to it. This is a some¬ 
what paradoxical situation since the general rule dictates 
that quantum walks turn into random walks when deco¬ 
hered [13 , and our starting point is a process whose dis¬ 
tribution already coincides with the classical one: when 
t is large enough 

(61) 

where p, = 2p — 1, = 4p(l —p), and the factor 2 is due 

to the odd/even alternating nature of n —see Fig[5l 
For homogeneous, time-independent coins as, e.g., the 
Hadamard coin, the introduction of random measure¬ 
ments on both position and chirality of the quantum par¬ 
ticle is enough to induce the shift from the original bal¬ 
listic behavior of the system, with a variance that scales 


as to the classical behavior of a random walk with a 
variance that grows linearly with time, and whose exact 
expression is a function of the measurement probability 
per unit of time, q [13. This approach does not suffice in 
the present case, due to the aforementioned existing link 
between chiral and positional degrees of freedom. 

Hence, we have to remove any residual correlation in 
the system to ensure the validity of the central limit the¬ 
orem: We will perform measurements of the location of 
the particle at random times and, after having pinned 
down its position, we will restore the chirality of the sys¬ 
tem to its initial state [13 . With this method, our QW 
turns into a sum of independent random variables which 
are not identically distributed, since the coin operator is 
still site- and time-dependent. Q 

This lack of homogeneity in the coin operator implies 
that the random walk obtained through decoherence will 
be a biased random walk, even if p = 0.5. This bias 
may lead to limiting probability density functions that 
are Gaussian, but it is well known [^ that depending 
on the strength of persistence, bimodal distributions may 
appear: Hence, in general, we will not recover the same 
statistical properties of the process whose evolution we 
were mimicking. 

We can observe all these phenomena in Fig. [S] Small 
measurement probabilities, q <C 0.1, result in normal dis¬ 
tributions with variance reduction. After that point ap¬ 
proximately, q > 0.1, variance starts increasing with g, 
and eventually becomes larger than the original value. 
This happens for measurement probabilities q > 0.25 
—see the collapse of the plots corresponding to g = 0 
and g = 0.25 in Fig. [SI Finally, for values of the measure¬ 
ment probability of about g > 0.7, bimodal distributions 
do appear. 


^ In some practical implementations of the discrete-time QW, the 
role of the time variable is assumed by an auxiliary spatial di¬ 
mension, what breaks in practice temporal homogeneity Is^l . 
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Figure 6. (Color online) Probability mass function p(n,t) for 
t = 100, p = 0.5, and different values of the measurement 
probability q. Data were obtained by averaging 1,000,000 dif¬ 
ferent realizations of the process, except for the case q = 0.00 
that corresponds to the coherent walk. Labels are ordered to 
indicate increasing variance. 


VII. CONCLUSION 

Inspired by the fact that quantum walks are univer¬ 
sal computation primitives, and thus they can solve any 
problem that can be tackled by a general-purpose com¬ 
puter, I looked for a particular instance that reproduced 


the statistical features of a random walk. 

The aim was to design a non-trivial version of the 
discrete-time quantum walk on the line with exactly the 
same probability of site occupation as the classical pro¬ 
cess at any time scale, not as a byproduct of the asymp¬ 
totic loss of coherence in the quantum evolution. Along 
the text, I have proved that one possible way to get the 
desired behavior is through the introduction of an inho¬ 
mogeneous, time-dependent coin operator. 

The correlation level shown by both components of the 
wave function is so high that one can use it to restore the 
system to the same state previous to a measurement of its 
chirality. This perfect reversion can be performed with 
the only aid of unitary operators whenever one knows the 
output of the measuring process. Moreover, the analysis 
of the entropy of entanglement between positional and 
chiral degrees of freedom shows that the latter informa¬ 
tion loses significance as time increases. 

Finally, this restoring procedure can be seen as a sim¬ 
ple protection mechanism against accidental degradation 
of the coherence of the quantum state, but it can lead to 
some other yet undiscovered interesting implications. 
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